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FORTRAN PROGRAM FOR PREDICTING OFF-DESIGN 
PERFORMANCE OF RADIAL-INFLOW TURBINES 
by Charles A. Wasserbauer and Arthur J. Glassman 
Lewis Research Center 

SUMMARY 

A FORTRAN IV program for calculating off-design performance of a radial-inflow 
turbine is presented. The program uses a one-dimensional solution of flow conditions 
through the turbine along the mean streamline. The loss model accounts for stator, 
rotor, incidence, and exit losses. Other program features include consideration of 
stator and rotor trailing-edge blockage and computation of performance to limiting 
loading. Overall turbine geometry and design-point values of efficiency, pressure 
ratio, and mass flow are needed as input information. The output includes perform- 
ance and velocity-diagram parameters for any number of given speeds over a range of 
turbine pressure ratio. Included in this report are the engineering equations, the pro- 
gram listing, and the input information and output listing for a sample turbine problem. 

The experimental performance of two radial- inflow turbines is compared with the 
results from the computer program presented herein and with those from a previously 
used program. The overall computed results from the program of this report show a 
marked improvement over those from the previously used program and good agreement 
with experimental data. 


INTRODUCTION 

A procedure for predicting the off- design performance of a radial-inflow turbine is 
described in reference 1. This procedure is based on a one -dimensional solution to 
determine the flow conditions through the turbine along the mean streamline. The 
FORTRAN IV computer program for this procedure is given in reference 2. This pro- 
gram uses turbine geometry and design-point performance as input and computes the 
performance and velocity- diagram parameters over a range of speed and pressure 
ratio. The results from this program, however, did not correlate as well as desired 



with the experimental results of some subsequently tested turbines. 

The basic loss-model assumptions used in this program were reexamined to see if 
alternative assumptions would yield a better correlation with experimental data. Ac- 
cordingly, the computer program described in reference 2 was modified. In addition, 
new features were added to the program to account for trailing-edge blockage and to 
compute the conditions for pressure ratios at or beyond stator and/or rotor choke. 

These features give more flexibility to the program and offer a greater range for off- 
design turbine analysis. 

This report presents the modified FORTRAN IV program which computes overall 
performance for a radial- inflow turbine. The revised loss assumptions and new pro- 
gram features are described. Also included in the report are the engineering equations, 
the program listing, and a complete description of input and output including a sample 
turbine problem. The performance of two radial- inflow turbines (refs. 3 and 4) was 
computed by using the program of this report and that of reference 2. The performance 
predicted by each program is compared with experimental data for the two turbines. 


METHOD OF ANALYSIS 


A one-dimensional solution of flow conditions through the turbine along the mean 
streamline is the basis for this analysis. Two independent variables are assumed for 
each calculated performance point. A value of rotor-inlet-tip speed Ug is chosen, 
and for each speed a range of stator-exit critical-velocity ratios (V/V cr ) is assumed. 


Symbols are defined in appendix A, and a typical turbine with station nomenclature is 
shown in figure 1. The equations used in the program are given in appendix B. The 
computer uses these equations in essentially the order listed. 

There are two modes of operation for this program - the design mode and the off- 
design performance mode. Before the off-design performance mode can be used, two 
loss determinants, the stator total-pressure ratio Pj/Pq and the rotor loss coefficient 
K, must be known. In the design mode the program automatically determines these 
loss determinants at the design point by means of a search routine. By repeated itera- 
tions the computer finds a pair of values for Pj/Pq and K which satisfy the input de- 
sign values of mass flow rate, efficiency, and pressure ratio. These values of Pj/Pq 
and K are then used in the off-design performance mode. If desired, the values of 
Pj/Pq and K can be determined from the user's own loss procedures and used in the 
off-design performance mode. The off-design performance mode computes all per- 
formance parameters over the entire range of rotor speeds and stator-exit velocity 
ratios requested. 
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Turbine Losses 


The results from the computer program of reference 2 did not correlate well with 
the experimental results from the turbine tests of references 3 to 6. The three main 
losses were (1) a stator loss, (2) an incidence loss, and (3) a rotor loss. It was de- 
cided to reexamine and modify the models used to represent these losses to see if a 
better correlation could be obtained. 

Stator loss . - In reference 1 the loss in kinetic energy across the stator is propor- 
tional to the average kinetic energy in the blade row and is represented by the equation 


y2 _ 

L _ V l,id V 1 

s — 


= K c 


VQ + Vj 

2gJ 


( 1 ) 


The stator loss coefficient Kg is determined from design-point performance and then 
is assumed to be constant for the off-design calculations. 

In order to check this assumption, stator loss coefficients were calculated from 
unpublished data obtained with the turbines of references 3 to 6. All four of these tur- 
bines had static-pressure taps located just inside the stator trailing edge. Therefore, 
stator performance could be calculated by using these static-pressure data along with 
the design stator-exit blade angle and the measured values of mass flow rate, stator- 
exit area, and total temperature. The calculated stator loss coefficients are given as 
a function of stator-exit critical- velocity ratio in figure 2. As shown, the stator loss 
coefficient decreases significantly with increasing velocity ratio over the range of data. 
Thus, the assumption of a constant stator loss coefficient, as made in reference 1, 
does not seem to be valid. 

Another way to express stator loss is as the ratio of stator-exit to turbine-inlet 
total pressure. Figure 3 shows the variation of stator-exit to turbine-inlet total- 
pressure ratio with stator-exit velocity ratio as obtained from the data of references 3 
to 6. The stator total- pressure ratio for each turbine was fairly constant over the 
range of velocity ratio tested. Total- pressure -ratio level ranged from about 0. 98 to 
0. 99 for most of the data. This examination and analysis indicated that a constant 
stator total-pressure ratio would be a better model for the stator loss, and thus it was 
incorporated into the program. 

Incidence loss . - Minimum incidence loss does not occur at zero incidence angle 
with respect to the rotor blade, but at some optimum flow angle <p (fig. 4). The cal- 
culation of the optimum flow angle follows the method of Stanitz (ref. 7). This method 
was used in reference 2 and is described in appendix B as equations (B8a) to (B8d). 

The flow angle at the rotor inlet is /3g, and the incidence angle is defined as 
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In reference 2 the incidence loss is assumed to be the component of relative velocity 
normal to the angle <p : 


L m = 


W 3 sin2i 3 
2gJ 


(3) 


In the program of this report, different variations in loss with positive and negative 
incidence were used, and equation (3) was changed to 


L IN = 


wjj sin n l 3 
2gJ 


(3a) 


A value of 2 was used for the exponent n for negative incidence. However, for positive 
incidence, a value of 3 for the exponent n gave a better correlation with experimental 
data. 

Rotor loss . - In reference 1 the viscous loss in the rotor is assumed to be propor- 
tional to the average kinetic energy in the blade row as calculated from the equation 


, wild -*4 
= 

K 2gJ 



(4) 


At high pressure ratios the level of rotor inlet velocity Wg seemed to have an exces- 
sive influence on the loss. Using the component of velocity in the direction of the opti- 
mum flow angle gave a better correlation with experimental data. In the present pro- 
gram, therefore, the rotor loss was calculated as 




Wl 


2gJ 



(4a) 


The value of K was approximately 0. 3 for all the turbines examined in this report. 
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Trailing-Edge Blockage 


In order to account for the effect of stator trailing-edge blockage, the analysis of 
reference 1 assumes a variable station where the flow from the stator is assumed to 
occupy the entire cylindrical flow area. Trailing-edge blockage at the rotor exit is not 
taken into account. 

The analysis of this report accounts for trailing-edge blockage at both the stator 
and rotor exits. The effect of blockage for each blade row is specified in terms of the 
ratio of the flow area just inside to that just outside the blade trailing edge. Figure 5 
shows the blade-row trailing-edge region, specifically the areas used in the blockage 
calculation. Both angular momentum and continuity are conserved when the conditions 
at these stations are calculated. 


Stator and Rotor Choke 


The program of reference 2 did not compute turbine performance at or beyond the 
stator and rotor choking points. The program of this report allows turbine performance 
to be computed at the stator and rotor choking points and at pressure ratios beyond 
choking to rotor blade limiting loading. The stator choke point is where (V/V cr ) = 1- 0, 

which is the point of maximum flow per unit area. For values of (V/V cr ) greater than 

1. 0, a new stator-exit flow angle ot^ is computed from the area required to pass the 
choking mass flow rate. 

In order to find the rotor choking point, an iteration is required to determine the 
value of (V/V cr )^ that maximizes flow per unit area at the rotor exit. Conditions up- 
stream of the rotor exit are then held fixed. As the velocity ratio (w/W cr ) is in- 
creased beyond the choking value, the exit flow angle (3^ is adjusted to pass the choking 
mass flow. The program is terminated at or close to blade limiting loading, where 

( w x //w cr) 4 =■ !• °- 


COMPARISON OF COMPUTED AND EXPERIMENTAL RESULTS 

This section compares the experimental results obtained with two radial-inflow 
turbines to the results obtained by the analytical procedures described in reference 2 
and this report. The two radial-inflow turbines are those of references 3 and 4. The 
results of the comparison are presented in terms of mass flow rate and total and static 
efficiency variations with pressure ratio and speed. 
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Mass Flow Rate 


Calculated and experimental variations of mass flow rate with turbine total- to 
static-pressure ratio for various speeds are compared in figure 6(a) for the reference 3 
turbine and in figure 6(b) for the reference 4 turbine. The data for the reference 4 tur- 
bine are unpublished air data, which were used herein because they covered a wider 
range of pressure ratios than the published argon data. For both turbines, the off- 
design values of mass flow rate computed by using the program of this report showed a 
significantly better correlation with the experimental data, especially at higher pressure 
ratios and lower speeds, than did the values computed by the program of reference 2. 
The poorest correlation between the computed and experimental values was in the region 
near choked flow (fig. 6(b)). In this region, the maximum deviation between computed 
and experimental mass flow rates was reduced from about 9 percent (program of ref. 2) 
to about 3^ percent by using the program of this report. 


Total Efficiency 

Calculated and experimental variations of total efficiency with blade-jet speed ratio 
for a number of speeds are compared in figure 7(a) for the reference 3 turbine and in 
figure 7(b) for the reference 4 turbine. In figure 7(a) the 90- and 110-percent speed 
lines are not shown in order to avoid overlapping of data. The off-design values of 
total efficiency computed by using the program of this report showed a significantly 
better correlation with the experimental data, especially in the case of the reference 4 
turbine, than did the values computed by the program of reference 2. As shown in fig- 
ure 7(b), the maximum deviation between computed and experimental total efficiencies 
was reduced from about 10 percent (program of ref. 2) to essentially zero by using the 
program of this report. 


Static Efficiency 

Calculated and experimental variations of static efficiency with blade-jet speed 
ratio for a number of speeds are compared in figure 8(a) for the reference 3 turbine and 
in figure 8(b) for the reference 4 turbine. Except for blade- jet speed ratios higher than 
the design value, the efficiencies computed by the program of reference 2 for all speeds 
were generally about 2 percentage points lower than experimental values for both tur- 
bines. The efficiencies computed by the program of this report over the same range of 
blade-jet speed ratio generally were within 1 percentage point of the experimental 
values. 
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The overall improvement in the correlation of calculated values of mass flow rate 
and efficiency with experimental data indicates that the turbine loss assumptions used 
in this report provide a better model than those used in reference 2. 


FORTRAN PROGRAM 
Program Input 

The program input consists of two title, or heading, cards followed by the required 
physical data and option indicators in NAMELIST format. The information contained in 
columns 1 to 60 of the title cards is printed as two lines of heading on the output listing. 
The two title cards, even if left blank, must be the first two cards in the data package. 
Two additional title cards must precede each different case being run in the same data 
package . 

The physical data and option indicators are input in data records having the 
NAMELIST name IN. All necessary physical data and option indicators must be inputted 
for the first case in the data package. For succeeding cases in the same data package, 
only those items being changed need be inputted. 

Options . - There are three sets of options that must be specified by the input. All 
three are specified by the variable MODE as described in the input variable list. The 
first option is the choice of units, SI or U. S. customary, to be used for input and output. 
The particular unit to be used for each variable is included in the input variable list. 

The second option is the choice of a mode of operation: either the design mode or 
the off-design performance mode. The design mode is used to automatically determine 
the stator total- pres sure ratio and rotor loss coefficient that yield design flow and effi- 
ciency at the design pressure ratio. The input variable ITEST, as described in the input 
variable list, is used to specify whether total efficiency and total- to total- pres sure 
ratio or static efficiency and total- to static-pressure ratio are used as the design 
values being matched. The off-design performance mode is used to compute the per- 
formance over the desired ranges of speed and pressure ratio. 

The third option provides for the choice of long or short output for the off-design 
performance mode. Long output is always given for the design mode. The long output 
includes complete velocity-diagram information in terms of critical- velocity ratios and 
angles, actual and equivalent overall performance parameters, and dimensionless de- 
sign parameters. The short output includes only certain of the equivalent performance 
and dimensionless design parameters. The exact output provided is described in the 
section Description of Output. 

I nput variables . - The input variables comprising NAMELIST IN are as follows: 
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AL1 


AO 

A1 

A3 

A5 

BL1 

BL4 

B4 

DELV 

DELY 

D2 

D3 

ETAD 

G 

ITEST 


MODE 


stator-exit blade angle, a j, deg 

2 2 

area at turbine inlet (scroll inlet), A Q , m ; ft 

area upstream of stator exit (circumferential area inside blade trailing edge), 
A v m 2 ; ft 2 

area upstream of rotor inlet (circumferential area outside blade leading edge), 
Ag, m 2 ; ft 2 

area downstream of rotor exit (annular area outside blade trailing edge), Ar, 
m 2 ; ft 2 

blockage factor at stator exit, ratio of area inside blade passage to area out- 
side blade passage (fig. 5) 

blockage factor at rotor exit, ratio of area inside blade passage to area out- 
side blade passage (fig. 5) 

rotor-exit blade angle, 0^, deg 

incremental value of stator-exit critical- velocity ratio (V/V cr )^ 

incremental value of rotor-exit critical- velocity ratio (W/W cr ) used after 
rotor choke 


stator-exit diameter, D 2 , cm; in. 

rotor-inlet diameter, Dg, cm; in. 

specified design value of static efficiency r] s or total efficiency 7] t (see 
ITEST) 

ratio of specific heat at constant pressure to specific heat at constant volume, 

r 

specifies which design values are used for PDSGN and ETAD: 

if ITEST=1, total- to static -pressure ratio Pq/p 5 and static efficiency 77 g 
are used 

if ITEST=2, total- to total- pres sure ratio Pq/p^ and total efficiency 77 1 
are used 

specifies which program option is used: 

MODE=0 yields the off-design performance mode in SI units and with long 
output 

MODE =1 yields the off-design performance mode in SI units and with short 
output 

MODE =2 yields the design mode in SI units 
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MODE =3 yields the off-design performance mode in U. S. 
and with long output 

MODE =4 yields the off-design performance mode in U. S. 
and with short output 

MODE =5 yields the design mode in U. S. customary units 

customary units 
customary units 

PD 

stator total- pressure ratio, Pj/Pg 


PDSGN 

specified design total- to total-pressure ratio pg/pg or total- to static- 
pressure ratio Pg/pg (see ITEST) 

P0 

n 

inlet total pressure, Pg, N/m ; psfa 


R 

gas constant, R, J/(kg)(K); ft-lb/(lb)(°R) 


TO 

inlet total temperature, Tq, K; °R 


U3 

rotor-inlet tip speed, Ug, m/sec; ft/sec 


U4U3 

ratio of rotor-exit mean speed to rotor-inlet tip speed, U4/U3 (Rotor-exit 
mean speed must correspond to rotor-exit mean velocity diagram. It is 
recommended that this be at the area-mean radius if no better value is 
available. ) 

VMAX 

final value of stator-exit critical-velocity ratio (v/V cr ) 


VI 

initial value of stator-exit critical-velocity ratio (V/V cr ) 


WD 

design value of mass flow rate, w, kg/sec; lb/sec 


XK 

rotor loss coefficient, K 


YMAX 

final value of rotor-exit critical-velocity ratio (w/W cr ) 


ZZ 

number of blades at rotor inlet 



Sample input . - Input sheets with the data used in computing the performance for 
the reference 3 turbine (figs. 6(a), 7(a), and 8(a)) are shown in tables I and n. Selected 
output obtained with this input is presented and described in the next section. Table I 
is for the design mode. Table n is for the off-design performance mode for speeds of 
100, 110, 90, 70, 50, and 30 percent of design. Each line of the input form shown in 
tables I and II represents one data card. The first two cards are the mandatory title 
cards, which can contain any desired message. The next four cards contain the turbine 
physical data and option indicators. The design-point quantities ETAD, PDSGN, and 
WD are included for the design mode (table I). The loss coefficients XK and PD deter- 
mined by the design-mode calculation are included in the off-design-mode (table n) data. 
Additional data sets may follow the first data set, as they do in table n, and need only 


9 


include those values that differ from previous case data. As shown in table n, only the 
speed is changed for the next five cases. 


Description of Output 

The design-mode output for the input of table I is shown in table HI. The top line of 
output is a program identification title that is automatically printed. The next two lines 
are the title card messages. The following three lines are the input variable values. 

The symbolism used to identify the output values is defined in terms of the engineering 
symbols in the list at the end of this section. Printed next are the computed values of 
stator total-pressure ratio and rotor loss coefficient followed by the design values of 
mass flow rate, efficiency, and pressure ratio. 

The remainder of the output in table in is divided into two parts. In the section 
VELOCITIES AND ANGLES, all the absolute and relative critical- velocity ratios and 
angles at various stations throughout the turbine are listed for the design stator-exit 
critical- velocity ratio. The section OVERALL PERFORMANCE gives all the perform- 
ance parameters computed by the program. 

Off-design performance mode output is shown in tables IV and V. Table IV shows 
the long output and presents the first three points obtained by using the input of table H. 
Table V shows the short output that would be obtained from the first data set of table H 
if MODE-1. The first six lines in both tables are the same as for the design-mode 
output. The remainder of the output in table IV gives the velocities, angles, and overall 
performance in the same format as previously described for the design-mode output. 

In table V the output is limited to certain of the overall performance parameters useful 
for defining the overall performance map. 

A list of the variable names as used for the output and their corresponding engi- 
neering symbols is given in the following table : 
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Variable 

Engineering 

Variable 

Engineering 

Variable 

name 

symbol 

name 

symbol 

name 

ALPHA- 1 

“i 

ETA-T 

\ 

VCR)3 

ALP HA- 2 

a 2 

GAMMA 

7 

V/VCR )3 

ALPHA-3 

“3 

1-3 

i 3 

VR/VCR )3 

ALPHA- 5 

Gfc 

K 

K 

VU/VCR)3 

AO 

D 

A o 

N 

N 

Al 

A i 

NS 

N s 

VCR)5 

A3 

A 3 

N/T 

N /V^o 

V/VCR)5 

A5 

A 5 

NU 

V 

VU/VCR)5 

BETA- 3 

j8o 

P 

p 

VX/VCR)5 

BETA- 4 

O 

04 

PD 

( p » /p °) des 

BETA- 5 

05 

P0, 

p 0 

w 

BL1 

Bl 

P0, /P5 

p o/ p 5 

W. F. 

BL4 

B 4 

P0./P5, 

p o/ p 5 

WN/DEL 

DEL-H 

Ah’ 

R 

R 

WT/P 

DEL-H/T 

Ah’/TJj 

TOR 

r 

10 6 r/p{,; 


DESIGN ETA-S 

^s, des 

TOR/P 

W/WCR)3 

DESIGN ETA-T 

^t, des 
( p 'o /p 5) H 


i44r/p^ 


DESIGN P0,/P5 

TO, 

\J 

T b 

WR/WCR)3 

DESIGN P0,/P5, 

W des e 

U3 

U 3 

WU/WCR)3 

DESIGN WT-FLOW 

U4U3 

U 4 / U 3 

W/WCR)4 

D3 

°3 

V/VCR)0 

^ V/V er) 0 

WU/WCR)4 

E Q— DE L— H 
EQ-N 

Ah’ 

eq 
N n 

V/VCR)1 

(vAcr), 

WX/WCR)4 

EQ-PO, /P5, 

(>w_ 

VR/VCR)1 

( V r/ V cr) 1 

W/WCR)5 

EQ-PO, /P5 

(■¥< 

eq 

VU/VCR)1 

( V u/ V cr) 1 

WX/WCR)5 

EQ-TOR 

V/VCR)2 

(v/v cr ) 2 

EQ-W 

ETA-S 

eq 

^s 

VR/VCR)2 

(V V cr ) 2 

ZZ 


Engineering 

symbol 


V cr, 3 
( V / V cr) 3 

( V r/ V cr) 3 

(V V crl 

V cr,5 

( V / V cr) 5 

(V V cr) 5 

(V V cr) 5 

W 

WF 

WN/6 

10 OOOW-y/Tj/pf); 
144W -^/Tq/Pq 
( W / W cr) 3 
( W r/ W cr) 3 
( W u/ W cr) 3 
( W / W cr) 4 
( W u/ W cr) 4 
( W x/ W cr) 4 
( W / W cr) 5 

( W x/ W cr) 5 

ZZ 


Messages 

If there is no solution for a particular turbine case in the design mode, a message 
will be printed - "NO SOLUTION FOR THIS CASE - CHECK Al, BL1, AL1 OR B4, 
BL4. " This will happen when the program cannot select values for the loss determi- 
nants Pj/Pg and K that wil1 satisfy the design input values of mass flow, efficiency, 
and pressure ratio. 

In the off-design performance mode for conditions beyond choking to blade limiting 
loading, a message will be printed after the last computed performance point - "LAST 
CASE IS APPROXIMATE LIMITING LOADING CASE. " 


11 



MAIN PROGRAM 


FORTRAN Variables 

The FORTRAN variables used in the main program are defined in the following 
table in terms of the engineering symbols, where available, or by descriptive termi- 
nology: 


Variable 

name 

Engineering 

symbol 

Variable 

name 

AB 

<P 

DELHOT 

A LX 

a i 

DELL 

ALl 



AL2 

«2 

DELTA 

AL3 

“3 

DELV 

AL5 

“5 


AX 

temporary storage 


AO 

A 0 

DELVl 

A1 

A 1 


A3 

A 3 

DELY 

A5 

A 5 


BL1 

Bl 

DHIDS 

BL4 

B 4 

DHIDT 

BX 

temporary storage 

DHTCR 

B3 

^3 

D2 

B4 

04 

D3 

B4X 

04 

ELAM 

B5 

05 

EPS 

CAPQ 

w/p 5 

EQPRS 

CASE 

variable for title 



message 

EQPRT 

COSAL1 

cos 

EQTOR 

COSAL2 

COS Ofg 

ET 

COSAL3 

cos c* 3 

ETAD 

COSB3 

cos 0 3 

ETAS 

COSB4 

cos 0 4 

ETAT 

CX 

temporary storage 

EX 

DELH 

1 Ah’ 

F 


Engineering 

symbol 

Ah'/Tg 

incremental value 

of ( W / W cr) 4 

6 

incremental value 

of (V/v„_\ 


Variable 

name 


of (V/Vcr)! 

incremental value 

of ( V / V cr) x 
incremental value 
of ( w / w cr ) 4 

Ah id,s 

Ah id,t 


Te/6 

temporary storage 
^s, des or ^t, des 


temporary storage 
integration variable 


Engineering 

symbol 


integration variable 
sinig 

integration variable 

r 

temporary storage 
g 

7+ 1 

7 " 1 

(7 ~ l)/(7 + 1) 
7/(7 + 1) 

7/(7 - i) 

(V'Vl 


(W^orl 


HI temporary storage 

I index variable 

IND integer variable controlling 

logical sequence in 
CONTEST 

IT index variable 

ITEST integer variable controlling 

proper input in GETK 
and SEEKPR 

J index variable 

K index variable 

L index variable 

M index variable 


12 



■■■■■■■ , «^“ l 1 


I I Hill II I I 


I 1 1 ■■ 


II 1 1 I II II II I 


II I II 


I I I I 


Variable 

name 

Engineering 

symbol 

MODE 

integer variable controlling 


all program options 

NAME 

variable for title message 

NCOUNT 

index variable 

P 

P 

PD 

p’i/po 

PDSGN 

("i^des " 

PR 

temporary storage (pres- 


sure ratio) 

PRS 

( p i/ p 5 ) 

PVCR 

( p ' v cr) 1 

PVK 

(P/P') cr 

PV1 

(pv/p-v cr ) i 

PV2 

(PV/P'V cr) 2 

PWC3 

(pw/p-w cr ) 3 

PWC4 

(P" W cr) 4 

PWVC3 

(P”w cr / P 'v cr ) 3 

PW4 

(p w/p-w cr ) 4 

PW4PW3 

(P" W cr , 4 /P" W cr,3) 

PO 

Po 

P3P3 

P 3 /P 3 

P4P3ID 

(P4’/P3) id 

P4P4ID 

(P4/P4’). d 

P5 

P5 

P5P0 

P 5 /P 0 

P5P0G 

(P5/Pb) (r ' 1)/y 

P5P0P 

P 5 /P 0 

P5P0PG 

(P 5 /p{,) (r ’ 1)/y 

P5P5 

(P’/P") 5 

P5P5P 

(P/P') 5 

Q 

temporary storage 


Variable 

name 

Engineering 

symbol 

QX 

(pv/p'v cr ) o 

Q1 

( T y T ”) 1/2 

R 

R 

RH05 

p 5 

RH05P 

p 5 

RH05RT 

(P/P% 

SCAPQ 

(Wi/p 5 ) 1/2 

SINAL1 

sin aj 

SINAL2 

sin a 2 

SINAL3 

sin 

SINB3 

sin 0 3 

SINB4 

sin j3 4 

SINB5 

sin j3 5 

SS 

temporary stora 

STHETA 

(0 )!/2 
v cr 7 

T 

0. 01745329 

THETCR 

0 cr 

TOL 

0. 00001 

TOP 

r / Po 

TOR 

r 

TO 

T o 

T3T3 

(T"/T’) 3 

T4T3 

T^’/T” 

T4T4 

x ,/ T m 

U3 

U 3 

U4U3 

V U 3 . 

VC5VC0 

^cr, 5 /^cr, 0^ 

VCR 

V cr,0 

VCR5 

^cr, 5 

VK 

1.0 

VMAX 

maximum value 


(WJ, 

VRVC1 

(V V cr), 


Variable 

name 

Engineering 

symbol 

VRVC2 

( V r /V cr) 2 

VRVC3 

(W) 3 

VU3T 

V u, 3, opt 

VU4VC3 

-l ^ 

CO 

VUVC1 

VUVC2 

(V V cr), 

VUVC3 

(V V cr) 3 

VUVC5 

(V V cr) 5 

VUWC4 

(V W cr) 4 

VUWC5 

(V W cr) 5 

VX 

temporary storage 
of variable VI VC 

VXVC5 

(V V cr) 5 

VOVCO 

( V/V cr)„ 

VI 

initial value of 

(V^cr), 

V1VC1 

i^cr). 

V2VC2 

( V / V cr) 2 

V3TU3 

^lopt 

V3VC3 

( V / V cr) 3 

V5VC5 

( V/V cr) 5 

V5WC5 

( V / W cr) 5 

WC4WC3 

W cr, 4 / W cr,3 

WCR3 


WCVC3 

( W c/ V cr) 3 

WD 

w des 

WK 

cr 

WNODE 

wNe/6 

WRWC3 

( W r/ W cr) 3 
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Variable 

Engineering 

name 

symbol 

WTHODE 

Me cr ) 1/2 e/5 

WTOP 

w(T 0 ) 1 / 2 /P 0 

WU3T 

w 

w u , 3, opt 

WUVC3 

( W u/ V cr) 3 

WUWC3 

( W u/ W cr) 3 

WUWC4 

( W u/ W cr) 4 

WUWC5 

( W u/ W cr) 5 

WXWC4 

( W x/ W cr) 4 

WXWC5 

( W x/ W cr) 5 

W1 

w 

W3VC3 

( w /v cr ) 3 

W3WC3 

( W / W cr) 3 

W3WC4 

W 3/ W cr,4 

W4WC4 

( W/W cr) 4 

W5WC5 

(w/w cr ) 5 

X 

iteration variable 


for velocity 


Variable 

Engineering 

name 

symbol 

XDEL 

index variable on 
estimated ve- 
locity in 
CONTIN 

XI3 

h 

XJ 

j 

XK 

K 

XN 

N 

XNOT 

N/(TM 1/2 

XNOTH 

N/e 1 / 2 
/ cr 

XNS 

N s 

XNU 

V 

XI 

iteration variable 
for velocity 

Y 

integration variable 
for velocity 

YCAL 

temporary storage 

YGIV 

temporary storage 

YMAX 

< w/w A™ x 

Z 

W.I J 5 

ZZ 

ZZ 

Z1 

Po/P5 


Program Listing 


A FORTRAN IV COMPUTER PROGRAM TO PREOICT DESIGN AND OFF 
DESIGN PERFORMANCE OF A RADIAL INFLOW TURBINE 

INTERNATIONAL SYSTEM OF UNITS 

M0DE=0 » IMPLIES OFF-DESIGN PERFORMANCE MODE WITH FULL OUTPUT. 

M0DE=1 t IMPLIES OFF-DESIGN PERFORMANCE MODE WITH SHORT OUTPUT. 

M0DE=2, IMPLIES DESIGN MODE TO DETERMINE THE VALUE OF K AND PD. 

U. S. CUSTOMARY UNITS 

MODE=3» IMPLIES OFF-DESIGN PERFORMANCE MODE WITH FULL OUTPUT. 

M0DE=4, IMPLIES OFF-DESIGN PERFORMANCE MODE WITH SHORT OUTPUT. 

MODE=5, IMPLIES DESIGN MODE TO DETERMINE THE VALUE OF K AND PD. 

DIMENSION NAME ( 10 > .CASE ( 10» 

NAMEL IST/IN/ G » PO , TO , AO , A 1 , A3 , A5 ,D3 , AL I , B4 , BLI ,BL4 ,PDS GN , I TEST , 
1J3.U4U3.R .Vl.DELV.VMAX.XK. WD.ZZ.MODF.PD. ET AD .DELY . YM AX .02 
SUBR(X,G)=( l.-( G- I.) /(G+l. >*X*X)** (l./(G-l.) >*X 
1 READ! 5» 100) ( NAME (I> , I =1 , 10) 

READ! 5 . 100) <CASE!I>. 1=1.10) 

REAOI 5 ? IN ) 

WR ITE( 6.200 )NAME. CASE 

WRITE! 6»400>G» P0»A1.AL1»B4»PD»R»T0»A3»BL1»BL4»XK»D3»A0»A5»U3» 
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1U4U3,ZZ 

IFCMODE.EQ. 1* OP # MODE .E Q. 4 > WR ITF<6.666> 

D EL V 1 = DEL V 
ALX=AL 1 
64X=94 
T=. 01745329 
TOL = .OT331 
AL 1 = T* AL 1 
B4=T*B4 

C0SAL1=C0S(AL1 ) 

S I NAL 1=S I N ( AL 1 ) 

C0SB4=C0S< B4) 

S IMB4=SIN( 84> 

X J = 778.029 

GR = 32.1741 

IT= l 

J= 1 

K “ 1 

L = 1 

M=1 

N= 1 

vx=o. 

XDEL= • 1 
G 1=G+ 1 • 

G 2 = G- 1 • 

G3=G2 / G 1 
G4= G/Gl 
G5 = G/G2 

BEGIN STATOR ANALYSIS 

IFCMODE.EQ. 3. OR .MODE .£Q.4.0R»M0DE.EQ.5> GO TO 2 
VCR=SQRT(2.*G4*R*T0> 

PVCR= P0*$QRT( 2.*G4/R/T0> 

GO TO 3 

2 VCR- SQRTC 2.*G4*GR*R*TQ> 

P VCR= PO* SQRTC 2.*G4*GR/R /TO) 

3 V1VC1=V1 
9 CONTINUE 

IFCV1VC1.GT.1. > GO TO 15 
NC0UNT=0 

I F ( MODE.EQ. 0 .OR.M ODE .EQ. l.OR. M0DE.EQ.3»0R. MODE .EQ« 4 ) GO TQ 10 
PD = 1.0 

4 K- 1 

SS=WD/PVCR/A1/C0SAL1/PD 
V 1 VC 1 = * 5 

5 X = V1VC1 

F = SUBR(X,G) -SS 

FP= f F+SS>/X-2**X*X/G1*C 1 • -G3* X*X 1 ** ( 1. /G2-l*> 

VlVCl^ X-F/FP 

IF(ABS<(V1VC1-X)/V1VC1).LT.T0L) GO TO 10 
GO TO 5 
15 VK = 1. 

PVK=SUBR< VK.G) 

WK=COS( ALX*T>*PVK *A1*PVCR*PD 
PV1=SUBR( V1VC1,G> 

CO SAL 1=WK/ 1 A 1*PD* PV1*P VCR ) 

AL 1=AC0S( COSAL 1 ) 

S I NAL 1 = SI N ( AL 1 ) 

GO TO 13 
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10 PV1=SU8R( VlVCltG) 

W 1=PV 1*A1*PVCR*C0 SAL 1* PD 
V0VC0=V1VC 1/6. 

QX=PV1*A1/A0*C0SAL1*PD 

12 X= VOVCO 

F= SUBR ( X, G ) -OX 

FP = ( F+QX )/X— 2.*X*X/G1*( 1.-G3* X*X ) ** ( 1./G2-1.) 
VOVCO= X-F/FP 

IF(ABS( ( VOVCO-X)/V1VC1).LT.TOL) GOTO 13 
GO TO 12 

13 VUVC1=V1VC1*SINAL1 
VRVC 1=V1VC 1*C0SAL 1 
VUVC2=VUVC 1 

STATOR EXIT CONDITIONS 

Q=PV1*C0SAL1 *BLl 
X*V1VC 1/3. 

61 F»I l.-G3*( X*X+VUVC2**2 > >** ( 1. /G2)*X-Q 
F1=(F+Q>/X 

FP=-2.*X*X/G1*F1**{ 1./G2-1. >+Fl 
X 1=X-F/FP 

IF<ABSC<X1— X)/X1) — TOL) 62.62.63 
63 X=X1 

GO TO 61 

62 VRVC2—X 

V2VC2=SQRT(VRVC2**2+VUVC2**2> 

PV2=SUBR(V2VC2»G> 

C0SAL2=VR VC2/V2VC 2 
AL2=AC0S<C0SAL2> 

SINAL2=SINtAL2) 

FREE STREAM SPACE CONDITIONS 

VUVC3=VUVC2*D2/D3 
Q* PV2*C0SAL2*D2/D3 
X*VUVC3/5. 

22 F = ( l.-G3*<X*X+VUVC3**2> >**U. /G2>*X-Q 
F1=(F+Q)/X 

FP=-2.*X*X/G1*F1**(1./G2-1.»+F1 
X 1=X-F /FP 

I F ( ABS I (X 1— X ) / X 1) —TOL ) 20,20,21 
21 X=X1 

GO TO 22 

20 V3VC3=SQRT (VUVC3* *2+X*X) 

C0SAL3=X/V3VC3 
AL3=AC0S(C0SAL 3) 

S INAL3=SIN( AL3) 

VRVC3— X 

ROTOR INLET CONDITIONS 

T3T3= ( 1.0— G3*C 2.0*U3*VUVC3/VC R— C U3 /VCR >**2>> 
WCVC3= SORT ( T3T 3) 

WUVC3 = VUVC3— U3/VC R 
WUWC3=WUVC3/WC VC3 
W3VC3=SQRT(WUVC3**2+VRVC3**2> 

W3WC 3=W3VC3/WC VC3 
WRWC3=SQRT C W3WC3**2— WUWC 3**2 ) 

P3P3=T3T3**(G5> 

PMVC3=T3T3**<Gl/2./G2> 

SINB3=WUWC3/H3WC3 
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B3=ASIN( S INB3) 

C0SB3=C0S C B3 ) 

V3TU3=1.-1.98/ZZ 
VU3T=V3TU3*U3 
WU3T= VU3T— U3 

AB=ATAN( ( WU3T/VCR >/VP,VC3) 

X I3=B3-AB 
FS =SIN<XI3) 

FC =COS( XI 3) 

WCP 3= WC VC 3* VCR 

T4T3=1.-G3*(U3/WCR3)**2* t l.-U4U3**2> 

HC4WC3=SQRT(T4T3» 

W3 WC4=W3WC 3/WC 4WC 3 
PW4PW3=T4T3**< U/G3/2. > 

PWC3=SUBR( W3WC3.G > 

YGIV=PWC3*<A3/A5/BL4 )*<COS B3/C0SB4) /PW4PW3 
IF(M0DE«EQ.2.OR.MODE*EQ.5> XK = .l 
C ROTOR ANALYSIS 

C 

Y= .3 

31 I ND= 1 

32 NCOUNT=NCOUNT+l 
CX=W3WC4**2*( XK*FC**2+FS**2) 

IF(B3.GT.AB) CX=W3WC4**2*( XK*FC**2+FS**3) 

A X= XK*Y*Y+CX 

R X = l.-G3*Y*Y 

EX= < i.-G3*<AX/BX >)**G5 

YCAL=SUBR<Y»G)*EX 

IF( IND.GE. 6.AND.ABSI ( YC4 L- YG I V > /YC AL ) . LE.TOL > GO TO 83 
CALL CONTIN< Y. YCAL.INO.M. YGIV.XOEL) 

I F ( IND.L T . 10) GO TO 32 
C IND=10 INDICATES CHOKEO FLOW 

I F { IND.EQ* 10) GO TO 33 

C IND = 11 INDICATES NO SOLU T I ON IS FOUND IN 130 ITERATIONS 

I F ( IND.EQ. 11) WR I TE ( 6» 2040 > 

I F ( NCOUNT • GE • 1000 ) WRITE (6.2050) 

33 V1VC1= V1VC1-DELV 

34 DELV= DELV/2. 

V1VC1= V1VC 1+DELV 
N = 2 

GO TO 9 

35 I T=2 

GO TO 82 

80 IF(M0DE.EQ.2.0R.N0DE.E0.5) GO TO 82 

I F < ABS( (V1VC1-VX) /V1VC1) .LT.TQL) GO TO 35 
VX=V1VC1 

IF(N.E0.2) GO TO 34 
82 W4WC4= Y 

PW4= SUBR ( Y * G ) 

P4P4ID=YGIV/PW4 

PWC4 =PVCR*PD*PWVC3*PW4PW3*P4P4ID 
85 CONTINUE 

WUWC4=W4WC4*S INB4 
W X WC4=W4WC4*C0SB4 
HA = WXWC4 

IF(HA.GT. 1. ) GO TO 56 

Q=PW4*C0SB4*BL4 

X=W4WC4/3. 
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70 F=( i*-G3*( X*X+HUWC4**2) /G2>*X-Q 
Fl=(F+Q)/X 

FP=—2« *X/G 1*F 1**( 2.-G>*X+Fl 
X1*X-F/FP 

I F ( AB S ( (X 1— X ) /X 1) — TOL) 71,71,72 
72 X = X1 

GO TO 70 

ROTOR EXIT CONDITIONS 


71 MXWC5=X 

W5WC5=X*X+WUWC4** 2 
W5WC5 = SQRT < W5WC 5) 

S I NB5=WUWC4/W5WC5 

VUWC4= 1. / WCVC3*U4U3*U3/VCR/WC 4WC3+ WUWC4 
VUWC5= VUWC4 
V5WC5=VUWC5**2+X*X 
V5WC5 = SQRT t V5HC 5) 

T4T4=1.+G3*(VUWC4**2-WUWC4**2I 

Q1*SQRT(T4T4) 

VXVC5=WXWC 5/Q1 
VUVC 5=VUWC5/Q1 
V5VC5*V5WC 5/Q1 
AL5=l»/T*ASINt VUVC5/V5VC5 ) 

B5=1./T*ASIN(SINB5> 

WUWC5=WUWC4 

VU4VC3=WUWC4*WC4WC3*WCVC3+U4U3*U3/VCR 
T 5T0- l.-2.*G3/VCR*(U3*VUVC3-U4U3*U3*VU4VC3> 

T5= T 0*T5T0 

VCR5 = VCR*SQRT(T5T0) 

P5P5 = T4T4**(G5) 

P4P3I0 = T4T3**<G5> 

P5P0 = PD*P3P3*P4P3ID*P4P4ID*P5P5 
P5=P0*P5P0 

P5P5P = <1.-G3*CV5VC5**2> )**(G5) 

P5P0P = P5P0*P5P5P 
VC5VC0 = < VCR5/VCR»**2 
P5P0G = P5P0**(G2/G> 

P5P0PG = P5P0P**( G2/G) 

Z = 1./P5P0 
Z1 = 1./P5P0P 
H1=1.-VC5VC0 
ETAS=H1/C 1.-P5P0PG) 

ETAT*Hl/( 1.-P5P0G) 

C AUTOMATIC DETERMINATION OF XK ANO PD 

C 

IF(M0DE.EQ.0.0R.M0DE.EQ.1.0R.MODE.EQ.3.OR.M0DE.EQ.4) GO TO 40 
ET= ETAS 

IFI ITEST.EQ.2) ET= ETAT 
IF(K.EQ.5> GO TO 30 
CALL GETMXK.ET.ETAD.tO 
IF( (ETAO-ET).GT.. 1) GO TO 90 
GO TO 31 
30 PR=Z1 

IFIITEST.EQ.2) PR=Z 
IFIL.EQ.5) GO TO 40 
CALL SEEKPRIPDSGN.PR.PD.L) 

GO TO 4 
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EXTRA OUTPUT CALCULATIONS 


40 CONTINUE 

IF(M0DE.EQ.3.OR.MODE.EQ.4.OR.MODE.EQ.5> GO TO 41 
THETR=G4*R*T0/48 247.36 
STHETA = SQRT(THETCR) 

DELTA =P0/ 101325. 

XNOTH= 100.*U3/3. 14159/D3/STHETA 

XN= XNOTH*STHETA 

DHTCR=H1*48247.36/G3 

DELH=DHTCR*THETCR 

DH IDS = ( l.-P5P0PG)*G5*R*T0 

DHIDT = C 1.-P5P0G >*G5*R*T0 

EQPRS = 1. /( 1.— DH IDS / 2 894 84. 2/THETCR)**3.5 
EQPRT = l./(l.-DHIDT/289484.2/THETCR)**3.5 
XNU=U3/SQRT<2.*DHIDS) 

WTOP = W1*100.**2*SQRT( TO) /PO 

T 0R= DELH*W1*9. 549274/XN 

TOP= TOR*1.E*06/P0 

ELAM = DELH/U3**2 

RH05P = P5P0*P0/VCR5**2*2.*G4 

P = DELH*W1/1000. 

GO TO 43 

41 THETCR=G4*R*T0/16141.4357 
STHETA = SQRT(THETCR) 

DELTA=P0/2116,22 
XN0TH=720.*U3/ 3. 14159/03/ STHETA 
XN= X NOTH* STHETA 

DHTCR=H 1*16141. 43 57/G3/XJ 

DELH=DHTCR*THETCR 

DH IDS = ( l.-P5P0PG)*G5*R*T0/X J 

DHIDT = ( 1.-P5P0G )*G5*R*T0/XJ 

EQPRS = 1. /C l.-DH IDS/124. 4808/THETCR) **3. 5 

EQPRT = l./Il.-DHIDT /1 24. 4808 /THETCR >**3.5 

XNU=U3/SQRT(2.*GR*XJ*DHIDS> 

WTOP=W 1*SQRT t TO )* 144. /PO 

TOR= DELH*W1/L. 12 164E-05/XN 

T 0P= T0R*144./P0 

ELAM = DEL H*GR*XJ /U3**2 

RH05P = P5P0*P0*2.*GR*G4/VCR5**2 

P = DELH*Wl*XJ/550. 

43 P R S= 1. /EQPRS 

EPS=< Gl/2. )**( G5> *.7395945/G 
WTH0DE=W1* STHETA* EPS /DELTA 
E QTOR = TOR*EPS/ DELTA 
WNODE = WTHODE*XNOT H 
DELHOT=DELH/TO 
XNOT= XN/SQRTITO) 

RH05RT = P5P5P**( l./G) 

RH05 = RH05RT+RH0 5P 
CAPQ = W1/RH05 
SCAPQ = SQRT(CAPQ) 

IF(MODE.EQ.3.0R.MODE.EQ.4.0R.MODE.EQ.5> GO TO 42 
XNS=XN*. 10472* SCAPQ/ (DHIDT**. 75) 

GO TO 45 

42 XNS=XN*SCAPQ/( ( X J*DHIDT) **. 75 ) 

45 CONTINUE 
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COMPLETE WRITE OUT 

IF<IT,EQ.3> GO TO 44 
AL 1=AL 1/T 
AL2=AL 2/T 
X 13 * XI3/T 
AL3=AL3/T 
B3*B3/T 
B4=B4/T 
44 CONTINUE 

IF(M0DE#EQ,1,0R,M0DE,EQ.4> WRITE(6,505» VI VC I » XNS» PRS , XNU, ETAT , 
1ETAS. OHTCR . WTHODE .EQPRS, EQPRT 
IF{M0DE.EQ.0.0R,M0DE,EQ.3> GO TO 46 
IF(MODE.EQ,1,OR.MODE.EQ,4) GO TO 48 
WRITEC6.509) PD.XK.WD 

I F C ITEST.EQ.l) WRITEC6,525) E T AD» PDSGN 
I F ( ITEST.EQ.2) WRITE(6,526> ETAO. PDSGN 
46 CONTINUE 

WRITE ( 6»500 )V1 VC 1 ,AL1 , V3VC3, W4WC4 , W5WC 5, VUVC 1, AL2, VUVC3, WUWC4, 
IWXWC5. VRVC1.AL3«VRVC3.WXWC4.V5VC5.V0VC0.B3.W3WC3.B4,VUVC5,V2VC2. 
2XI3.WUWC3.B5.VXVC5»VRVC2»VCR,WRWC3,VCR5,AL5 
WR ITE I 6» 50 1 ) WT MODE, EQPRS, Z l,WTOPtWl, OHTCR » ETAS ,XNU .DELHOT, DELH, 
IXNOTH. ETAT. XNS. XNOT* ELAM. EQTOR. EQPRT, Z. TOP, TOR tXN.P.WNOOE. T 5. P5 
I F ( MOD E ,E Q, 2 • OR • M ODE ■ E Q« 5 > GO TO 51 
48 CONTINUE 

IFCWXWC4.EQ, I, ) GO TO 59 
IF( J,EQ,2> GO TO 57 
IFI IT.GE.2) GO TO 55 
AL 1*ALX*T 
B4*B4X*T 

V1VC1=V1VC 1+DELV 
IFCVIVCI.GT.VMAX) GO TO 51 
GO TO 9 
51 AL1=ALX 

DELV* DEL VI 
B4* B4X 
GO TO 1 

55 I T*3 

GIV*YGI V*C0S(B4X*T> 

Y* Y+DELY 

IF(Y.GT.YMAX) GO TO 51 
53 W4WC4* Y 

PW4* SUBRCY.G) 

AX* XK*Y*Y+CX 
BX= 1,-G3*Y*Y 

P4P4ID* ( 1,-G3*(AX/BX) )**G5 
C0SB4*GIV/(PW4*P4P4ID> 

B4* ACOS( C0SB4) 

B4= -B4 
C0SB4* COS ( B4) 

S I NB4= SIN< B4) 

B4= B4/T 

IF(J.EQ*2) GO TO 58 
GO TO 85 

56 IF(J.EQ.2> GO TO 57 
Y = Y-DELY 

DELL* DELY/5. 

J= 2 

57 Y* Y+DELL 
GO TO 53 
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58 HB* W4WC4*C0SB4 
IF(ABSC(HB-HA)/HA I.LT..02) GD TO 59 
DELL* DELL/2. 

GO TO 85 

59 WR ITE ( 6*600) 

GO TO 51 

90 WRITE (6.300) 

GO TO 1 
STOP 

FORMAT STATEMENTS 


ERRDR MESSAGE NUMBER 1 
100 FORMAT ( 10A6) 

200 FORMAT C 1H1 , 35X,43HNASA RADIAL INFLOW TURBINE COMPUTER PROGRAM/ 
124X. 1 0A6/24X. 10A6 ) 

300 FORMAT ( /10X, 54HN0 SOLUTION FOR THIS CASE - CHECK A1,BL1,AL1 OR B4, 
1BL4) 

400 FORMAT ( 7HK GAMMA. 2X.G13. 5.7H PO. .2X.G13.5 .6H A1.2X.G13.5. 

12X»7HALPHA-1,2X,G13.5,7H BETA-4, 2X ,G13. 5 ,5H PD,2X,G13.5/ 


27H 

R , 2X, G13. 5, 7H 

TO, , 2X,G13.5.6H 

A3.2X.G13.5.2X. 

37H 

BLl.2X.G13. 5. 7H 

BL4.2X.G13.5.5H 

K. 2X ,G 13.5 / 

47H 

D3, 2X,G13. 5, 7H 

AO , 2 X, G13# 5 , 6H 

A5 ,2X, G13. 5, 2X , 

57H 

U3.2X.G13. 5.7H 

U4/U3.2X.G13.5.5H 

ZZ.2X.G13.5 ) 


500 FORMAT ( / 50X , 21HVEL0C I TI E S AND ANGLES/ 

112H V/VCRU »G13.5» IX » 1 1HA L PHA-1 ,G1 3. 5, IX , 1 1HV/ VCR) 3 , 

AG13.5. 1X.11HW/WCR )4 . G 13. 5 . IX. 1 1HW/WC R) 5 .G13.5/ 

212H VU/VCR)1 ,G13.5, IX , 11HALPHA-2 ,G13. 5, IX ,1 1HVU/VCR) 3 , 

AG13.5, IX, 11HWU/WCR)4 , G13. 5 , IX, 1 1HWX/WCR )5 ,G13.5/ 

31 2H VR /VCR ) 1 .G13.5. IX. 11HALPHA— 3 .G13. 5. IX . I 1HVR/VCR ) 3 . 

AG 13. 5, IX, 1 1HWX/WCR)4 ,G13.5 , IX, 1 1HV/VCR) 5 ,G13.5/ 

412H V/ VCR > 0 ,G13.5, IX. UHB ETA-3 . G1 3. 5. IX • 1 1HW/WCR) 3 , 

AG13.5. IX. 11HBETA-4 »G 13. 5 • IX. 1 1HVU/VCR )5 ,G13.5/ 

512H V/VCR) 2 , G 13. 5 , 1X,11HI— 3 ,G13. 5, IX ,1 1HWU/WCR ) 3 , 

AG13.5.1X.11HBETA-5 .G13. 5. IX. 1 1HVX/VCR )5 .G13.5/ 

61 2H VR /VCR ) 2 ,G 13. 5 , IX , 1 1HVCR)3 , G1 3. 5, IX , 1 1HWR/WCR ) 3 , 

AG13.5, 1X»11HVCR>5 , G13. 5 , IX, 1 1HALPH A-5 ,G13.5) 

501 FORMAT ( 50X • 19H0VERALL PERFORMANCE/ 

112H EQ-W ,G13»5, IX, 11HE Q— PO, /P5 , G1 3. 5 , IX , 1 1HP0 , / P5 

BG13.5, 1X,11HWT/P , G 13. 5 , IX , 1 1HW ,G13.5/ 

212H EQ-DEL-H .G 13. 5. IX. 11HETA-S .613.5. IX .11HNU 

BG13.5, IX, 11HDEL-H/T , G13. 5 , IX , 1 1HDEL-H ,G13.5/ 

312H EQ-N .G13.5. IX.11HETA-T .G13. 5. 1X.1 1HNS 

BG13.5. 1X.11HN/T .G13. 5 , IX, 1 1HW.F. ,G13.5/ 

412H EQ-TOR ,G 13. 5, IX , 1 1HE Q-PO , /P5 , , G13. 5, IX ,1 1HP0 , /P5 , 

BG13.5. 1X.11HT0R/P .G 13.5 . IX. 1 1HT0R .G13.5/ 

512H N ,G 13. 5, IX , 1 1HP , G13. 5, IX , 1 1HWN/DEL 

BG13.5, 1X,11HT5, , G13. 5 , IX , 1 1HP5 , ,G13.5) 

505 FORM AT (10G 13.5) 

509 FORMAT ( 13X , 2HP0, G 15. 6/ 14X , 1HK , G15. 6/1X , 14HDE SI GN WT-FL0W.G15.6) 

525 FORMAT ( 3X, 12HDE SIGN ETA-S ,G1 5.6/2 X, 13HDES IGN P0./P5.G15.6) 

526 FORMAT (3X. 12HDESI GN ETA-T .G1 5.6/ 1 X. 14HDES IGN PO ,/ P5 , , G15. 6 ) 

600 FORMAT ( /10X,41HLAST CASE IS APPROXIMATE LIMITING LOADING) 

666 FORMAT (7HK V1VC1. 8X.2HNS • 1 IX. 9HEQ— P5/P0. .4 X.2HNU. 1 1X.4HET AT .9X • 
14HETAS.9X. 5HDHTCR ,8X,6HWTH0DE ,7X ,9HEQ-P0 , / P5 ,4X, 10HEQ-P0 , / P5, ) 
2040 FORMAT! 10X.44HNO SOLUTION COULD BE FOUND IN 100 ITERATIONS) 

2050 FORMAT! 10X . 69HI TERATI ON PROCEDURE HAD TO BE RESTARTED TO AVOID NEG 
1ATIVE TEMPERATURE/15X, 67HRESTART PROCEDURE WAS ABORTED AFTER 1000 
2T0TAL NUMBER OF ITERATIONS) 

END 


Subroutine GETK (XK, ETA, ETAD, K) 


Subroutine GETK varies the value of the loss coefficient K by false positioning 
until the design value of efficiency is met. 


XK rotor loss coefficient, K 

ETA computed value of efficiency from main program 
ETAD design value of 77 ^ or 77 g 
K indicator used in method of false positioning 


SUBROUTINE GETK ( X K. ETA, ETAD. K > 

GO TO (100, 101), K 

100 XK1= XK 

D I F 1= ETA-ETA0 
K = 2 

X K= XK + .005 
R E T URN 

101 XK2= XK 

0IF2 = ETA-ETAD 
I F { 0 I F2*D I F 1 ) 104,103,102 

102 X K 1= X K 2 
XK= XK+.C05 
D I F 1= 0 1 F 2 
RE T UPN 

103 XK = XK2 
K = 5 
RETURN 

1C4 XK= -C XK2-XK1)/(D IF2— D IF 1 ) *DI F 1+XK 1 
K = 5 

RETURN 

END 


Subroutine SEEKPR (PDSGN, PR, PD, L) 

Subroutine SEEKPR varies the value of the stator total-pressure ratio Pj/Pq by 
false positioning until all design specifications are met. 

PDSGN design value of Pg/Prj or Pg/p 5 

PR computed value of turbine pressure ratio from main program 

PD value of stator total- pressure ratio Pj/Pg 

L indicator used in method of false positioning 
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SUBROUTINE SEEK PR ( PDSGN.PR.PD, L) 

GO TO ( 10 , 1 1 ) , L 

10 PDI=PD 

D I F 1= PR-PDSGN 
L = 2 

P 0=P0- • 005 
RETURN 

11 PD2=PD 

DIF2 = PR-PDSGN 
I F ( D I F2*D I F 1 ) 14,13,12 

12 P D1 = P D2 
PD=PD-.005 
D I F 1= DIF2 
RETURN 

13 PD=PD2 
L = 5 
RETURN 

14 PD= -( PD2-PD1) /<D IF2-DIF 1)*DIF 1+PD1 
L = 5 

RETURN 

END 


Subroutine CONTIN (XEST, YCALC, IND, JZ, YGIV, XDEL) 


Subroutine CONTIN is a curve-fitting routine which is described in detail in refer- 
ence 8. It is used to determine the rotor-exit velocity ratio value needed to satisfy 
continuity at the rotor exit. 

XEST value of estimated velocity Y 

YCALC mass flow parameter based on estimated velocity 

IND index to control next iteration in CONTIN and to indicate when a choked-flow 

solution has been found 

JZ index variable 

YGIV input mass flow parameter 

XDEL maximum permitted change in estimated velocity Y per iteration 


SUBROUTINE CONTIN < XEST, YCALC. I ND , J Z. YG I V .XDE L> 

C 

C — CONTIN CALCULATES AN ESTIMATE OF THE RELATIVE FLOW VELOCITY 
C — FOR USE IN THE VELOCITY GRADIENT EQUATION 
C 1 

DIMENSION X( 3) • Y( 3) 

NC ALL = NCALL+1 

IF ( IND.NE. 1.AND.NCALL.GT.100) GO TO 160 
GO TO <10*30.40.50.60,110.150) .IND 
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C— FIRST CALL 
10 MCALL = 1 
XORIG * XEST 

IF IYCALC.GT.YGIV.AND. JZ.EQ.l) GO TO 20 
IND = 2 
YC 1) = YCALC 
XI 1> * 0. 

XEST = XEST+XOEL 
RETURN 
20 IND * 3 

Y<3) * YCALC 
XI 3) * 0. 

XEST * XEST-XDEL 
RETURN 

C— SECOND CALL 
30 IND = 4 

Y<2) » YCALC 
XI 2) » XEST— XORIG 
XEST = XEST+XDEL 
RETURN 
40 IND = 5 

Y< 2) * YCALC 
XI 2) = XEST-XORIG 
XEST * XEST-XDEL 
RETURN 

C — THIRD OR LATER CALL - FIND SUBSONIC OR SUPERSONIC SOLUTION 
50 YI3) * YCALC 

XI 3» = XEST-XORIG 
GO TO 70 
60 YC1) = YCALC 

XI 1» * XEST-XORIG 

70 IF IYGIV.LT.AMIN1IYI1),Y<2>,Y(3))> GO TO (120, 130), JZ 
80 IND = 6 

CALL PABC I X»Y» APA »BPB,CPC ) 

DISCR * BPB**2— 4.*APA*ICPC— YGI V) 

IF IDISCR.LT. 0*> GO TO 140 

IF I ABS(400.*APA* ICPC-YG IV) ) • LE.BPB**2 ) GO TO 90 
XEST * — BPB-SIGNI SQRTCDISCR) , APA) 

IF I JZ.EQ. l.AND.A PA.GT.O.. AND. Y(3).GT.Y(1) ) XEST = -BPB+ 
1SQRTI DISCR) 

IF < JZ.EQ.2.AND.APA.LT.0. > XEST * -BPB-SQRT< DI SCR) 

XEST = XEST/2./APA 
GO TO 100 

90 IF < JZ.EQ.2.AN0.BPB.GT.0. > GO TO 130 
ACB2 * APA/BPB*ICPC-YGIV)/BPB 
IF <ABS(ACB2).LE. l.E-8) ACB2=0. 

XEST * -CCPC-YG IV )/BPB*< 1.+ACB2+2* *ACB2**2) 

100 IF (XEST.GT.XI3)) GO TO 130 
IF (XEST.LT.XI 1 ) ) GO TO 120 
XEST * XEST+XORIG 
RETURN 

C — FOURTH OR LATER CALL - NOT CHOKED 
110 IFIXEST-XORIG.GT.Xf 3)) GO TO 130 
I FI XEST— XORIG.LT. XI 1) ) GO TO 120 
YC2) * YCALC 
XI 2) = XEST-XORIG 
GO TO 70 

C — THIRD OR LATER CALL - SOLUTION EXISTS, 

C — BUT RIGHT OR LEFT SHIFT REQUIRED 
120 IND * 5 
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C — L EFT SHIFT 

X E ST = X( 1 J-XDEL+XOR IG 
X 0 SHF T = XEST-XORIG 
XORIG = X E ST 
VI 3) = Y (2) 

X(3> = X(2>-XOSHFT 
Y ( 2) = Y ( 1 ) 

X< 2) = X( l J-XOSHF T 

RETURN 
130 IND = 4 

C--RIGHT SHIFT 

X E ST = X{ 3 l+XDEL+XOR IG 
XOSHFT = X EST-XOR IG 
XDRIG = XE ST 
Y ( 1 > = Y(2) 

X<1> = X(2>-XOSHFT 

Y ( 2) = Y< 3 ) 

X(2> = X(3)-X0SHFT 
RETURN 

C--THIRD OR LATEP CALL - APPEARS to BE CHOKED 
140 XE ST = -BPB/2./APA 
IND = 7 

IF (XEST.LT.XI 1»> GO TO 120 
IF(XEST.GT.X(3) ) GO TO 130 
X E ST = XEST + XORIG 
RETURN 

C — FOURTH OR LATER CALL - PROBABLY CHOKED 
150 IF ( YCALC.GE.YGI V ) GO TO 110 
IND = 10 
RETURN 

C — NO SOLUTION FOUND IN 100 ITERATIONS 
160 IND = 11 
RETURN 
END 


Subroutine PABC (X, Y, A, B, C) 

Subroutine PABC calculates the coefficients A, B, and C of the parabola 

o 

y = Ax + Bx + C passing through three given X, Y points supplied by subroutine 
CONTIN. 



noon 


SUBROUTINE PA8C(X ,Y,A,B,C > 

— PABC CALCULATES COEFFICIENTS A.B.C OF the PARABOLA 
— y=A*X**2+B*X+C» PASSING THROUGH The giv c n x,y points 

DIMENSION X(3),Y(3) 

Cl = X ( 3 ) — X ( 1) 

C 2 = ( Y(2)-Y( 1) )/ t X ( 2>-X( 1) > 

A = (C1*C2-Y(3>+Y (1) >/Cl/(X(2 >-Xl3) I 

B = C 2- C X ( 1)+X( 2) )*A 

C * Y< 1>-Xl 1>*B-X ( 1>**2*A 

RETURN 

END 


Lewis Research Center, 

National Aeronautics and Space Administration, 
Cleveland, Ohio, June 17, 1975, 

505-04. 
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APPENDIX A 


SYMBOLS 

A area, m^; ft^ 

B blockage factor 

Cj dimensional constant, 200; 720 /tt 

Cg dimensional constant, 1; 360/rr 

Cg dimensional constant, 1000; 550 

D diameter, cm; in. 

9 

g conversion constant, 1; 32. 1741 ft/sec 

Ah' specific turbine work, J/kg; Btu/lb 

AhJ^ ideal turbine work based on inlet-total- to exit- static-pressure ratio, J/kg; 
Btu/lb 

Ah?d j- ideal turbine work based on inlet- total- to exit-total- pressure ratio, J/kg; 
Btu/lb 

i incidence angle, deg 

J mechanical equivalent of heat, 1; 778. 029 ft-lb/Btu 

K rotor loss coefficient, dimensionless 

Kg stator loss coefficient, dimensionless (ref. 1) 

L kinetic energy loss, J/kg; Btu/lb 

N turbine speed, rad/sec; rpm 

N g specific speed, dimensionless; rpm (ft®/^)/sec^^ 

n incidence loss exponent 

P power, kW; hp 

p pressure, N/m ; psfa 

R gas constant, J/(kg)(K); (ft-lb)/(lb)(°R) 

T absolute temperature, K; °R 

U blade speed, m/sec; ft/sec 

V absolute velocity of gas, m/sec; ft/sec 

W gas velocity relative to rotor, m/sec; ft/sec 
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WF work factor, eq. (B38) 
w mass flow rate, kg/ sec; lb/sec 
ZZ number of rotor blades at rotor inlet 

a absolute gas angle (angle between absolute velocity vector and meridional plane, 
positive when tangential velocity component is in direction of wheel velocity), 
deg 

/3 relative gas angle (angle between velocity vector relative to wheel and meridional 
plane, same sign convention as for a), deg 

r torque, N-m; in. -lb 

y ratio of specific heat at constant pressure to specific heat at constant volume 

6 ratio of turbine inlet total pressure to U. S. standard atmospheric pressure, 

eq. (B22) 

e function of y used in relating parameters to those using air inlet conditions at 
U. S. standard sea- level conditions, eq. (B23) 

77 0 efficiency based on ratio of inlet total to exit static pressure 

efficiency based on ratio of inlet total to exit total pressure 

& cr squared ratio of critical velocity at turbine inlet to critical velocity at U. S. 
standard atmospheric temperature, eq. (B21) 

v blade-jet speed ratio, eq. (B37) 

q o 

p gas density, kg/m , lb/ft 
<p optimum rotor flow angle, deg 
Subscripts: 

cr condition corresponding to V/V = 1 
des design 

eq air equivalent (U. S. standard sea level) values 

id ideal 

IN incidence 

max maximum 

opt optimum 

R rotor 

r radial component 
28 
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S stator 

u tangential component, positive when in direction of wheel velocity 
x meridional component, component in plane containing axis of rotation 

0 station at turbine inlet 

1 station immediately upstream of stator exit 

2 station immediately downstream of stator exit 

3 station immediately upstream of rotor inlet 

4 station immediately upstream of rotor exit 

5 station immediately downstream of rotor exit 
Superscripts: 

’ absolute total state 

M total state relative to rotor 

* U. S. standard sea-level air conditions (temperature, 288. 15 K (518. 67° R); 

pressure, 101 325 N/m (2116.22 psfa), specific-heat ratio, 1.4; gas constant, 
287.04 J/(kg)(K) (53.35 (ft-lb)/(lb)(°R)) 
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APPENDIX B 


EQUATIONS USED IN PROGRAM 


The analytical procedure involves a step-by-step solution of the flow conditions 

through the turbine along a mean line. Thus, at the rotor exit, the flow conditions and 

velocity diagrams are those at the mean radius, which could be a flow mean or an area 

mean. The two independent variables that are fixed for any given calculation point are 

the rotor-inlet-tip speed ratio (u/V ) and the stator-exit critical-velocity ratio 

3 

(V/V ) . The analytical procedure and the equations used are outlined in the following 
paragraphs . 


Stator Analysis 

The total temperature is assumed to be constant for the first four stations. Thus, 

V cr, 0 = V cr, 1 = V cr, 2 = V cr, 3 
The mass flow per unit area is expressed as 


pV 

p'V 
F v cr 



y - 1 

y + 1 



(Bl) 


For an input value of Pj/Pq and the assumed value of (V/V ) , the continuity relation 
at station 0 is 



p i/pv \ 

PoWcr^ 


cos a. 


(B2) 


Equation (Bl) is substituted into equation (B2), which can then be solved iteratively for 

( v/v cr) 0 - 
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The mass flow rate is then computed as 


w 


p i/pv \ 

PoV'Vcry' 


( p,V cr) 0 A l 


COS Q!j 


(B3) 


where (p'V cr ) 


0 


is evaluated from the known inlet conditions of Pq and Tq by using 


P' = 


P' 

RT' 


and 


V 


cr 


i/-^- gRT' 

V 1 


For values of V/V cr greater than 1. 0, the choking value of mass flow rate is used to 
calculate a new stator-exit flow angle: 


cos a i = 


w 


cr 


p i/pv \ 

PoWcr^ 


(P’ V cr) 



(B4a) 


At station 1 the geometry of the velocity diagram gives 



(B4b) 


(B4c) 


and since 


Station 2 conditions are determined by assuming that (p'V cr ) = (p’V cr ) 

D 0 = D, , (V /V = (V /V ) . The continuity relation between stations 1 and 2 is 
z v u cry^ v u cr /2 


written as 
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[ pV ) B 1 cos a 1 =[ — ) cos a.cy 
1 \P7, 


cry 


cry 


where 


and 


B 1 = — 
A 2 


pV 




ipT 


cos a 2 = <1 - — 


cr> 


y + 1 


cry 


u 


cry 


i/(r-l) 




cry 


The geometry of the velocity diagram gives 


cr> 




cr/ 


cry 


1/2 


and 


<22 = cos 


(-) 

_1 \ V cr/ 2 


cry 


Equations (B5a) to (B5e) are solved iteratively to determine (V r /V cr ) , (V/V cr ) 

2 2 


a 2‘ 


The conditions at station 3 are determined by assuming that 


(P' V cr) 2 


(B5a) 


(B5b) 


(B5c) 


(B5d) 


(B5e) 


, and 
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and 



The continuity relation between stations 2 and 3 is given as 



D c 


cos 0*2 = 



COS Q!g 


where 


r 



Hi- 


r - 1 

y + 1 





From the geometry of the velocity diagrams, 



and 


a 


3 



Equations (B6a) to (B6e) are solved iteratively to determine (V r /V cr ) , (V/V cr ) 

3 3 


a, 


3* 


(B6a) 


(B6b) 


(B6c) 


(B6d) 


(B6e) 


, and 
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Rotor Analysis 


The relations between relative and absolute parameters at the rotor inlet are 
by the following four equations: 




11 

T' 


tt\(y+l)/2(y-l) 


The velocity- diagram geometry gives 



given 


(B7a) 


(B7b) 


(B7c) 


(B7d) 


(B7e) 


(B7f) 


(B7g) 
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I 


'W. 


u 


i3g = sin 


-1 


\W 


cr/ 


W 


IW 


cr/ 


The optimum rotor- inlet flow angle cp is calculated as follows: 


'u, 3, opt _ i _ !• 98 


U 0 


ZZ 


r , . = u, A u ’ 3,opt> ) 

u, 3, opt 3 \ U 3 / 


W u, 3, opt V u, 3, opt " ^3 


cp = tan 


(1 u,3, 

-A V cr, 




The rotor incidence angle is 


cr/ 


h ~ @3 ~ v 


Since the rotor mean radius decreases from inlet to exit, there is a relative total- 
temperature drop expressible by the following equation: 


rp f » 

A 3 


y 


— f— f 

+ 1 \ W cr/„ 


1-1 




This allows the following rotor-exit parameters to be calculated: 


(B7h) 

(B8a) 

(B8b) 

(B8c) 

(B8d) 

(B9) 

(BlOa) 
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(BlOb) 


n »? /ffttvy/ly - !) 

p 4,id _/ 1 4 \ 

p 3 \ T 3/ 


W 

W 


cr, 

cr, 


4 

3 



(BlOc) 


f p ' W cr) /T ,,x(y+l)/2(y-l) 

V cr/ 4,id _ j 4 \ 

(P" w cr) 3 V$) 


(BlOd) 


The rotor-exit conditions are calculated by using the continuity equation between sta- 
tions 3 and 4 


(pAW r ) 3 = (PAWJ^ 


which is written as 



( pW \ A 3 cos ^3 
^ MW cr/ 3 A 4 cos ^4 

nil (P*^) 

p 4 V cry 4, id 
P4,id (‘“" w cr). 


(Bll) 


where = B^Ag. 

Everything on the right side of equation (Bll) is known except P 4 /P 4 id , which is 
the relative total pressure recovery for the rotor. It can be expressed as follows: 




(B12) 
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Equation (Bl) is substituted, in terms of relative quantities, into equation (Bll). And 
the values of (\V/W cr ) and P 4 '/P 4 ld are determined by an iteration procedure with 


equations (Bll) and (B12). 
rotor exit are held fixed. 


After the rotor choke point, conditions upstream of the 
As the velocity ratio (W/W cr ) is increased beyond the 


choking value, the exit flow angle /3 4 is adjusted by the following equation: 


cos 0 4 = 



(B13a) 


The velocity- diagram geometry gives 



(B13b) 



(Bl3c) 


Station 5 conditions are determined by assuming that (p M W cr ) = (p" w cr )^ and 

fw /W ) = (W /W„_} . The continuity relation between stations 4 and 5 is given as 
\ u' cr/ a \ \i cr/p- 


f pw 

tP" W, 


B 4 cos jS 4 = 


cr/ 


/ PW \ 

V^’ W cr L 


cos 


% 


(B14a) 


where 


B 


4 



(B14b) 


and 
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l/(y-l) 


r 


/LpW_\ COS h 4 - ^ 

\p" w cr/L ] >- +1 


'M 2 /M 2 


^ W er / 5 V w cr / 4 


2J 


The geometry of the velocity diagram gives 



% 




(Bl4c) 


(B14d) 


(Bl4e) 


Equations (B14a) to (B14e) are solved iteratively to determine the values of (W x /W cr ^ , 
(W/W cr ) 5 , and 0,. 

The relations between absolute and relative parameters at the rotor exit are given 
by 


where 



- 1 


y - i 

y + 1 



y - 1 

y + 1 



(B15a) 



jy \ _/T , \ y// ^ y ~ 1 ^ 
P'V 5 Vt ”/ 5 


(B15b) 


(B15c) 
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and 



(B15d) 


With the assumption that (V u /W cr ) = (V u /W cr ) and (W x /W cr ) = (V x /W cr ) , the 

4 5 5 5 

geometry of the velocity diagram gives 



(B16a) 


(B16b) 


(B16c) 


(B16d) 


Overall Turbine Performance 


The turbine overall total-temperature ratio is given by 


T 5 

— = 1 


o /r - 1\ f U 3 V u, 3 U 4 V u, 4' 


y + 1 / 


V 


cr, 3 


cr, 3/ 


(B17a) 
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where 



and the critical velocity at the turbine exit is 



The overall turbine total- to total- pressure ratio is given by the equation 

MUM /P4,id\/ P4 \/p'\ 

Po \P0APV3\P3 /\P 4 ,id/'P ,, /5 

The total- to static-pressure ratio at the turbine exit is obtained from 

y/(y- 1) 




The overall turbine total- to static-pressure ratio is then 


P5 

PO 



The turbine total and static efficiencies are obtained from 


*t = 


1 





p .Ur-l)/r 


(B17b) 


(B17c) 


(B18) 


(B19a) 


(B19b) 


(B20a) 
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(B20b) 


1 - 


= 


i-p 

Vo/ 


,(y-l)/y 


The following equations define the additional performance parameters which appear 
in the output listing: 


~ y + 1 


— RT o 


cr 


(t^t rt T 


5 = 


eb 


»** 


e = 


= 0.739594 5 ^y + Vj /' 1 ) 


N 


C 1 U 3 


eq 


D 3< e cr> 1/2 


Ah' =/Lz_\rZo 4-11 

eq T 4 , 


w. 


eq 


w(0 cr ) 1/2 € 

5 


CgJw Ah' 
eq N 6 


(wN) =■ 2^ 
eq 5 


(B21) 

(B22) 

(B23) 

(B24) 

(B25) 

(B26) 

(B27) 

(B28) 
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Ah' o 
id, s 


.-(2 

w 


{y-\)/y 


- 1 J 


Ah; 


id, t 


/ p ^\(y-l)/y 

U/ 


RT n 

y U 

- 1 J 




_ fy Z l V j Ah id, t 


yRT * 


~l-y/(y- 1 ) 


cr 



Ah' = Ah' 6 nr% 
eq cr 

N = N (6 ) 1//2 

eq v cr' 


r 



6 

e 


u 3 

< 2 S J 4 h id , S > 1/2 

WF - Si Ah ' 

U 2 

u 3 

p _ Ah' wJ 
C 3 


(B29) 

(B30) 

(B31) 

(B32) 

(B33) 

(B34) 

(B35) 

(B36) 

(B37) 

(B38) 

(B39) 
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WT _ 

p = p6 “ 


DEL-H _ Ah* 
T ~T’ 


N N 


T ( T ^/ 2 


TOR _ r_ 
p "pj, 


/T£N 

t 5- t oL7 

Vo/ 


p - pb (S 


, _ p 5 


Pk = 


RT' 


^5 


-*fcl 


1/y 

5 


N s = 


%) 

(J 


1/2 


0. 75 


(B40) 

(B41) 

(B42) 

(B43) 

(B44) 

(B45) 

(B46) 

(B47) 

(B48) 
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TABLE I. - SAMPLE INPUT FOR DESIGN MODE 
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| FORTRAN STATEMENT i IDENTIFICATION 


5 6 7 S 9 10 U 12 13 lb 15 16 17 18 ' 19 20 21 22 23 26j25 26 27 28 29 3d 31 32 33 36 35 36|37 38 39 60 61 42 l 63 644S 46 47 6849 50 St 52 55 54 55 56 57 58 59 60 ftl 62 63 66 65 6667 68 69 70 71 72 73 76 75 76 77 78 79 80 1 


i , 12.|62-CM iPI A TURBINE ..... . . . , , , . . 

[ 10 0, PERCENT SPEED I 


,$IN_ PO = 172 3 6 8 .9, T 0 = 1 1 4 4 . , 44 , A Q = . 0120176, Al=. 003 8 008 , A3 = . 0 , 03 9 7 8 , 9 , ,, A , 5 = .00445 , 4 , 3 
A L 1 =72.47,B 4 = -5 6 . 86 , , ,U3,=|2 3 7 . 9 52 4,,,U 4,U3,= . 5 6 7 7 6 5 , R =j 9 , 9 ,.,l 976 , G = 1 .j666,7 , i , , t I , , 

BlT=| . 9546 8 BL 4 =.937 18 , ZiZ = 22 , , V 1, = . 3 0 , DELV = , . , 0 5 , , VMAft = , ■ , 8 , 0 ; , DELY =[ . 05 , . ,4 0,,, 
DjP = 12.6 2,3,8 , ,pj?^12 . 9 9 97 , , MQD R , = Qe,T ^ = . 9, 1,3, , RDSfiN = 1, . 74 , ITRjS T =2 , , ;wb,= .,3 ,3,87,9, , $ , J , , 
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TABLE n. - SAMPLE INPUT FOR OFF-DESIGN PERFORMANCE MODE 


FORTRAN STATEMENT 


7 8 9 10 II 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30( 31 32 33 34 35 36 37 38 39 40 41 42 43 46 45 46 47 48 49 50 51 S2 53 54 55 56 57 58 59 6oj 61 62 63 64 65 66167 68 69 70 71 72 


73 74 75 76 77 78 79 80 


JIN 

AL1 

t-- I- H 

BL1 


D 3 = 


$1 N 


$IN 


$IN 


12 . 


62 -CM 


DIA TURBINE 


100 


PERCE 


NT SPEED 


0 = 172 3 
2. 47, B 


95168 


6 8 .9 , 

4 =-56 , 
tB,L4 = .91 


0 = 11441 


01 


4_4, Ad=.. 
8_6_ a U3 =12 3 r. 9^2 4 J[ 
371 8, ZZ=2 2 . ,V1 = 


U4U[3 = . + 5 6 7| 
S OjDEL V, 


. , 6 2 3 8, , 


02 = 12 . 


12. 62-CM 


ID I A TURBINE 


110 PERCE 


NT SPEED 




3=261 . 7476 $ 


12 . 62-CM DIA TUR BI NE 


90 PERCENT SPEED 


3=214 .1572 $ 


, 1, 2, . 62 -CMf Djl A T UR B INE 
7 0 PERCENT SPEElT^ 


U3, = 16 6 . 56 67 $ 


12 .62-CM DIA 


TURBINE 


2 0 


28 6 9,9, , ( PD_=^9 8 7^83_ J_ 




= .003 8 


0 08 , A 3 = . 003 9 7 89 


99. 1 9 7 6 . G = 1 . 66 


7 6 5 j R=[ 

. J) f 5* VMA^ = . , 8, 0 . DftL Y , = . 05 


mm 


5 = 


. 0 


50 PERCENT 


SPEED 


$ 1 N 


U 3=118.97 62 $ 


12. 62-CM DIA 


TURBINE 


30 PERCENT 


SPEED 


$1 N 


U 3 =71.3 8 , 57 $ 


12 3 4 5 6 7 8 9 10 II 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 3 5 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 6B 69 70 71 72 1 


73 74 75 76 77 78 79 80 
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TABLE IE. - OUTPUT FOR DESIGN MODE 


NASA RAD l At INFLOW TURBINE COMPUTER PROGRAM 
12.62-CM DIA TURBINE 
100 PERCENT SPEED 


GAMMA 

1.66670 

PO, 

0.17 237E+06 

A 1 

0.38008E-02 

ALPHA- 1 

72.4700 

BETA-4 

-56.8600 

PD 

500000000000 

R 

99.1976 

TO, 

1144.44 

A3 

Q.39789E-02 

BL1 

0.95168 

BL4 

0.93718 

K 

500000000000 

D3 

12.6238 

AO 

0.12018E-01 

A5 

0.44543E-02 

U3 

237.952 

U4/U3 

0.56776 

ZZ 

22.0000 


PD 0.9873 36 

K 0.286840 

DESIGN WT-FLOW 0.338790 

DESIGN ETA-T 0.913000 

DESIGN P0./P5. 1.740000 


VELOCITIES AND ANGLES 


V/VCRH 

0.60480 

ALPHA- 1 

72.4700 

V/ VCR) 3 

0.62055 

W/WCR )4 

0.42154 

W/WCR) 5 

0.41351 

VU/VCR) 1 

0.57671 

ALPHA-2 

73.2895 

VU/VCR >3 

0.59389 

WU/WCR ) 4 

-0.35297 

WX/WCR) 5 

0.21541 

VR/ VCR ) 1 

0-18217 

ALPHA-3 

73.1421 

VR /VCR ) 3 

0.17996 

WX/WCR )4 

0.23045 

V/VCR ) 5 

0.22205 

V/VC R) 0 

0.49333E-0 1 

BETA-3 

-11.8554 

W/WCR)3 

0.19253 

BETA-4 

-56.8600 

VU/VCR) 5 

0.37859E-01 

V/VCRI2 

0.60214 

1-3 

5.67598 

WU/WCR )3 

-0. 39554E-01 

BETA-5 

-58.6058 

VX/VCR)5 

0.21860 

VR/ VCR) 2 

0.17314 

VCR )3 

376.707 

WR/WCR >3 0.18842 

OVERALL PERFORMANCE 

VCR) 5 

340.824 

ALPHA-5 

9. 81656 

EQ-W 

0.22000 

EQ-P0,/P5 

1.68890 

PO ,/P5 

1.79477 

WT/P 

0.66492 

W 

0.33879 

EQ-DEL-H 

35013.4 

ETA-S 

0.86975 

NU 

0.69152 

DEL-H/T 

44.9930 

DEL-H 

51491.8 

EQ-N 

494.764 

ETA-T 

0.91298 

NS 

0.97688E-02 

N/T 

17.7359 

W. F. 

0.90941 

eq-tor 

148.673 

EQ-PO, /P5, 

1.64445 

PO ,/P5 , 

1.73997 

TOR/P 

1610.76 

TOR 

27 7.644 

N 

599.999 

P 

17.4449 

WN/DEL 

108.849 

T 5* 

936.601 

P5, 

99064.4 


-a 


TABLE IV. - LONG OUTPUT FOR OFF-DESIGN PERFORMANCE MODE 


NASA RADIAL INFLOW TURBINE COMPUTER PROGRAM 
I2.62-CM DIA TURBINE 
100 PERCENT SPEED 


GAMMA 

1* 

66670 

PO. 0* 1 723 7E+06 

A1 0.38008E-02 ALPHA-1 

72. 4700 

BETA-4 

-56.8600 

PD 0.98784 

R 

99 

.1976 

T0» 1144 

• 44 

A3 0.39789E-02 BL 1 

0.95168 

BL4 

0.93718 

K 0.28684 

D3 

12 

• 6238 

AO 0. 1 2018E-01 

A5 0.44543E-02 U3 

237.952 

U4/U3 

0.56776 

ZZ 22.0000 






VELOCITIES 

AND ANGLES 





v/vcru 


0.30000 

ALPHA-1 

72.4700 

V/VCR) 3 

0.30765 

W/WCR ) 4 

0.18127 

W/WCR ) 5 

0. 17792 

VU/ VCR) 1 


0.28607 

ALPHA-2 

73.2733 

VU/VCR )3 

0.29459 

WU/WCR )4 

-0.15179 

WX/WCR) 5 

0.92B32E-01 

VR/VCR) 1 


0. 90362E— 0 1 

ALPHA-3 

73.2404 

VR/VCR )3 

0. 88714E-01 

WX/WCP ) 4 

0.99100E- 

-01 V/VCR >5 

0. 23649 

V/VCR) 0 


0.27291E-01 

BETA-3 

-75.2549 

W/WCR)3 

0.34739 

BETA-4 

-56.8600 

VU/VCP )5 

0.21762 

V/VCR )2 


0.29871 

1-3 

-42.6024 

WU/WCR >3 

-0.33595 

BETA-5 

-58.5503 

VX/VCR )5 

0. 92548E-01 

VR/ VCR) 2 


0.B5970E-01 

VCR) 3 

376.707 

WR/WCR >3 

0.88418E-01 

VCR ) 5 

366. 180 

ALPHA-5 

66.9618 






OVERALL PERFORMANCE 





EQ-W 


0.12178 

EQ-PO » /P5 

1.25303 

PO ,/P5 

1.27858 

WT/P 

0.36807 

W 

0. 18754 

EQ-DEL-H 


10634.9 

ETA-S 

0.58860 

NU 

1.03221 

DEL-H/T 

13.6661 

DEL-H 

15640.0 

EQ-N 


494.764 

ETA-T 

0.68226 

NS 

0. 12937E-01 

N/T 

17.7359 

W.F. 

0.27622 

EQ-TOR 


24. 9971 

EQ-PO, /P5, 

1.21377 

PO , / P5 , 

1.23436 

TOR/P 

270.823 

TOR 

46.6815 

N 


599.999 

P 

2.93309 

WN/DEL 

60.2536 

T 5. 

1081.37 

P5. 

0. 13964E+06 






VELOCITIES 

AND ANGLES 





V/VCRll 


0.35000 

ALPHA- 1 

72.4700 

V/VCR) 3 

0.35895 

W/WCR )4 

0.21500 

W/WCR ) 5 

0.21102 

VU/ VCR ) 1 


0.33375 

ALPHA-2 

73.2751 

VU /VCR 13 

0.34368 

WU/WCR >4 

-0. 1 8003 

WX/WCR ) 5 

0.11009 

VR/VCR) 1 


0.10542 

ALPHA-3 

73.2298 

VR/VCR >3 

0.10357 

WX/WCR )4 

0.11754 

V/VCR >5 

0.22218 

V/VCR )0 


0.31446E-01 

BETA-3 

-70.2194 

W/WCR)3 

0.30739 

BETA-4 

-56.8600 

VU/VCP) 5 

0. 19303 

V/VCR)2 


0.34849 

1-3 

-41.4568 

WU/WCR >3 

-0.28926 

BETA-5 

-58.5551 

VX/VCR ) 5 

0. 11002 

VR/VCR) 2 


0.10029 

VCR) 3 

376.707 

WR/WCR >3 

0.10403 

VCR ) 5 

362.253 

ALPHA-5 

60.3183 






OVERALL PERFORMANCE 





EQ-W 


0.14031 

EQ-PO , /P5 

1.29942 

PO ,/P5 

1.33130 

WT/P 

0.42407 

W 

0.21607 

EQ-DEL-H 


14524.7 

ETA-S 

0.69588 

NU 

0.96036 

DEL-H/T 

18.6645 

DEL-H 

21360.4 

EQ-N 


494.764 

ETA-T 

0.77581 

NS 

0.12228E-01 

N/T 

17.7359 

W.F. 

0.37725 

EQ-TOR 


39.3343 

EQ-PO, /P5, 

1.26365 

P0./P5. 

1.29061 

TOR /P 

426.155 

TOR 

73.4559 

N 


599.999 

P 

4.6153 8 

WN/DEL 

69.4213 

T 5, 

1058.30 

P 5, 

0. 13356E+ 06 






VELOCITIES 

AND ANGLES 





V/VCR ) 1 


0.40000 

ALPHA-i 

72.4700 

V/VCR) 3 

0.41026 

W/WCR ) 4 

0.25008 

W/WCR ) 5 

0.24544 

VU/ VCR ) 1 


0.38142 

ALPHA-2 

73.2772 

VU/VCR )3 

0.39278 

WU/WCR >4 

-0.20941 

WX/WCP )5 

0. 12802 

VR/VCR) l 


0.12048 

ALPHA-3 

73.2172 

VR/VCR) 3 

0.11846 

WX/WCR ) 4 

0.13672 

V/VCR ) 5 

0.21080 

V/VCR) 0 


0.35422E-01 

BETA-3 

-63.6238 

W/WCR13 

0.26994 

BETA-4 

-56.8600 

VU/VCR)5 

0. 16728 

V/VCR ) 2 


0.39827 

1-3 

-37.9869 

WU/WCR )3 

-0.24184 

BETA-5 

-58.5610 

VX/VCRJ5 

0. 12827 

VR/VCR ) 2 


0.11460 

VCRI3 

376.707 

WR/WCR )3 

0.11993 

VCR ) 5 

358.262 

ALPHA-5 

52.5177 






OVERALL PERFORMANCE 





EQ-W 


0.15804 

EQ-PO, /P5 

1.35205 

P0./P5 

1.39174 

WT/P 

0.47763 

W 

0.24337 

EQ-OEL-H 


18435.5 

ETA-S 

0.77129 

NU 

0.89744 

DEL-H/T 

23.6900 

DEL-H 

27111.8 

EQ-N 


494.764 

ETA-T 

0.83787 

NS 

0. 11633E-91 

N/T 

17.7359 

W.F. 

0.47883 

EQ-TOR 


56.2318 

EQ-PO. /P5. 

1.31074 

P0./P5. 

1.35341 

TOR/P 

609.226 

TOR 

105.012 

N 


599.999 

P 

6.59807 

WN/DEL 

78.1904 

T5, 

1035.11 

P5, 

0.12736E+06 



TABLE V. - SHORT OUTPUT FOR OFF-DESIGN PERFORMANCE MODE 


NASA RADIAL INFLOW TURBINE COMPUTER PROGRAM 
12.62-CM DIA TURBINE 
100 PERCENT SPEED 


GAMMA 

1.66670 

PO. 

0. 17237E+06 

A1 

0.38008E-02 

ALPHA-1 

72.4700 

BETA-4 

-56.8600 

PD 

0.98784 

R 

99.1976 

TO, 

1144.44 

A3 

O.39789E-02 

BL1 

0.95168 

BL4 

0.93718 

K 

0.28684 

03 

12.6238 

AO 

0. 1 2018E-01 

A5 

0.44543E-02 

U3 

237.952 

U4/U3 

0.56776 

ZZ 

22.0000 


V1VC1 

NS 

EQ-P5/P0, 

NU 

ETAT 

ETAS 

DHTCR 

wthode 

EQ- PO ,/P5 

EQ-P0,/P5 

0.30000 

0. 12937E- 01 

0.79806 

1.03221 

0.68226 

0.58860 

10634.9 

0.12178 

1.25303 

1.21377 

0.35000 

0.12228E-01 

0.76957 

0.96036 

0.77581 

0.69588 

14524.7 

0.14031 

1.29942 

1.26365 

0.40000 

0.11633E-01 

0.73962 

0.89744 

0.83787 

0.77129 

18435.5 

0.158Q4 

1.35205 

1.31874 

0.45000 

0.11106E-01 

0.70763 

0. 84072 

0.87750 

0.B2169 

22379.6 

0.17485 

1.41317 

L. 38098 

0.50000 

0. 10630E- 01 

0.67325 

0.78877 

0.90064 

0.85235 

26373.0 

0.19064 

1.48533 

1.45245 

0.55000 

0.10196E-01 

0.63616 

0.74065 

0.91145 

0.86735 

30438.1 

0.20532 

1.57192 

1.53575 

0.60000 

0.98037E- 02 

0.59608 

0.69565 

0.91311 

0.86996 

34606.9 

0.21878 

1.67762 

1.63416 

0.65000 

0. 94640 E- 02 

0.55348 

0.65396 

0.90984 

0.86434 

38907.1 

0.23093 

1.80673 

1.74969 

0.70000 

0.91646E-02 

0.50642 

0.61352 

0.90D78 

0.84947 

43444.3 

0.24169 

1.97464 

1.89244 

0.75000 

0.89057E- 02 

0.45030 

0.57117 

0.88408 

0.82184 

48495.8 

0.25097 

2.22075 

2.08615 

0. 79767 

0.87530E-02 

0.34555 

0.50395 

0.84101 

0.74373 

56375.0 

0.25838 

2.89392 

2.51398 

0.79767 

0. 87904E- 02 

0.32379 

0.49136 

0.82970 

0.72255 

57612.4 

0.25838 

3.08840 

2. 61 149 

0.79767 

0.88588E- 02 

0.30207 

0.47914 

0.81733 

0. 70049 

58739.5 

0.25838 

3.31047 

2.71498 

0.79767 

0.89605E-02 

0.28066 

0.46737 

0.80367 

0.67765 

59721.8 

0.25838 

3.56308 

2.82441 

0. 79767 

0.90976E- 02 

0.25990 

0.45619 

0.78839 

0.65409 

60504.8 

0.25838 

3.84761 

2.93913 

0.79767 

0.92695E- 02 

0.24037 

0.44584 

0.77092 

0.62985 

61000.7 

0.25838 

4.16021 

3.05727 

0.79767 

0.94658E-02 

0.22317 

0.43681 

0.75035 

0.60520 

61060.3 

0.25838 

4.48087 

3.17450 

0.79767 

0. 96420E- 02 

0.21093 

0.43043 

0.72494 

0.58139 

60411.1 

0.25838 

4.74084 

3.28122 

0.79767 

0.96773E-02 

0.20992 

0.42990 

0.69101 

0.56152 

58490.2 

0.25838 

4.76379 

3.35536 

0.79767 

0.96545E-02 

0.21189 

0.43093 

0.68254 

0.55811 

57857.0 

0.25838 

4.71935 

3.36 237 


LAST CASE IS APPROXIMATE LIMITING LOADING 




Exit- to inlet-total- 
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Figure 7. - Comparison of calculated and experimental total efficiencies. 
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